Acoustical-Mode-Driven Electron-Phonon Coupling in Transition-Metal Diborides 
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We show that the electron-phonon coupling, Aq^, in the transition- metal diborides NbB2 and 
TaB2 is dominated by the longitudinal acoustical (LA) mode, in contrast to the optical E2g 
mode dominated coupling in MgB2- Our ah initio results, described in terms of phonon disper- 
sion, linewidth, and Xqy along V — A, also show that (i) NbB2 and TaB2 have a relatively weak 
electron-phonon coupling, (ii) the E2g linewidth is an order of magnitude larger in MgB2 than in 
NbB2 or TaB2, (iii) the E2g frequency in NbB2 and TaB2 is considerably higher than in MgB2, 
and (iv) the LA frequency at A for TaB2 is almost half of that of MgB2 or NbB2. 



The discovery of superconductivity in MgB2 ^ has 
renewed interests |, |, |, |, g, |, |, in find- 
ing superconductivity in similar materials such as simple 
metal diborides (BeB2, BeB2.7 5 Ipjl ), transition-metal di- 
borides {NbB2 I, I, TaB2 g %MoB2 §, ZrB2 §), 
borocarbides {LiBC H, |^), and other alloys jl^. Sur- 
prisingly, many of the materials mentioned above have 
not shown any superconductivity, those who are found 
to superconduct do so at a relatively low Tc < 10 K, 
and very likely, require hole doping and/or external pres- 
sure to show superconductivity. For example, the recent 
work of Yamamoto et al. [Q showing superconductivity 
in hole-doped NbxB2 is a case in point. Thus, of all the 
diborides, MgB2 seems to be in a class by itself with a 
Tc = 39 Kwhile all the other diborides have Tc < 10 K. 

Recent work on MgB2 P, [l||6|, |l|, ^ 
have unambiguously shown a very strong and anisotropic 
electron-phonon coupling in this system. In particular, it 
is found that electron/hole states on the cyHndrical sheets 
of the Fermi surface along T — A couple very strongly to 
the in-plane B — B bond stretching E2g phonon mode 
H, m, H, m, m. The E2g phonon mode coupHng 
gives rise to partial electron-phonon coupling constant 
Aqjy of the order of 2 — 3 Such a large partial Aqi, 

along r — A results in a relatively high superconduct- 
ing transition temperature of 39 K. Thus, the observed 
differences in the superconducting properties of MgB2 
vis-a-vis other diborides, specially the transition-metal 
diborides, can be better understood by comparing the 
electron-phonon coupHng along T — A in these systems. 
Taking the cue, we have studied from first principles (i) 
the phonon dispersion uj^^, (ii) the phonon linewidth jqi, 
and (iii) the partial electron-phonon coupling con- 
stant Xqi, along T — A in MgB2, NbB2 and TaB2 in 
PQ/mmm crystal structure. 

We have calculated the electronic structure of 
MgB2, NbB2 and TaB2 in P6/mmm crystal structure. 
We used optimized lattice constants a and c for MgB2 
and NbB2 ||2^, and experimental lattice constants (a — 
5.826 a. M., c = 6.130 a. w.) for TaB2 which are close to 
our optimized values of a = 5.792 a.u. and c = 6.143 a.u.. 
The lattice constants were optimized using the ABINIT 
program |S based on pseudopotentials and plane waves. 



For studying the electron-phonon interaction we used the 
full-potential Hnear response program of Savrasov ||2^, |2^ 
to calculate the dynamical matrices and the Hopfield pa- 
rameter, which were then used to calculate the phonon 
dispersion ujqi/, the phonon linewidth 7qy, and the par- 
tial electron-phonon coupling constant Xqi, along T — A 
in MgB2, NbB2 and TaBj. 

Based on our calculations, described below, we find 
that (i) in contrast to a strong and i?2g-mode domi- 
nated electron-phonon coupling in MgB2, the transition- 
metal diborides NbB2 and TaB2 have a relatively weak 
electron-phonon coupling which is dominated by the lon- 
gitudinal acoustical {LA) mode, (ii) the E2g phonon 
linewidth is an order of magnitude larger in MgB2 than 
in NbB2 or TaB2, and (iii) the i?2g phonon frequency 
in NbB2 as well as TaB2 is considerably higher than in 
MgB2 while the LA phonon frequency at A for TaB2 is 
almost half of that of or NbB2. 

Before describing our results in detail, we provide some 
of the computational details. As indicated above, the 
structural relaxation was carried out by the molecular 
dynamics program ABINIT ||2^ with Broyden-Fletcher- 
Goldfarb-Shanno minimization technique using Troullier- 
Martins pseudopotential for MgB2 and Hartwigsen- 
Goedecker-Hutter pseudopotentials ||2^ for NbB2 and 
TaB2, 512 Monkhorst-Pack ^ k-points and Teter pa- 
rameterization for exchange-correlation. The kinetic en- 
ergy cutoff for the plane waves was 110 Ry for MgB2 and 
140 i??/ for NbB2 and TaB2- The charge self-consistent 
full-potential LMTO calculations for electronic struc- 
ture were carried out with the generalized gradient ap- 
proximation for exchange-correlation of Perdew et al 
and 484 k-points in the irreducible wedge of the Brillouin 
zone. In these calculations, we used s, p, d and / orbitals 
at the Mg and Ta sites, and s, p and d orbitals at the 
Nb and B sites. The 2p state of Mg as well as the 5s and 
bp states of Ta were treated as semi-core states. In all 
cases the potential and the wave function were expanded 
up to Imax = 6. The muffin-tin radii for Mg, B, Nb and 
Ta were taken to be 2.4, 1.66, 2.3 and 2.5 atomic units, 
respectively. 

The calculation of dynamical matrices and the Hopfield 
parameters along T — A were carried out for 4 equidis- 
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Figure 1: The phonon dispersion ojq^ of MgB2 (upper panel), 
NbB2 (middle panel) and TaB2 (lower panel) for longitudinal 
acoustical (solid circles connected with full line) and optical 
E2g (open circles connected with dashed line) modes along F — 
A, calculated using the full-potential linear response method 
as described in the text. The lines connecting the points are 
only a guide to the eye. 



Table L The phonon frequencies (in meV) at F and A for 
the E2g and the LA modes for MgB2 calculated using the 
linear response method as described in the text as well as 
from previous work. 
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tant q-points for MgB2 and 7 equidistant q-points for 
NbB2 and TaB2. For Brillouin zone integrations we used 
a 12 X 12 X 12 grid while the Fermi surface was sampled 
more accurately with a 36 x 36 x 36 grid of k-points using 
the double grid technique as outlined in Ref. We 
checked the convergence of the relevant quantities by car- 
rying out Brillouin zone integrations using a 16 x 16 x 16 
grid of k-points with Fermi surface sampling done over 
48 X 48 X 48 grid. 

In Fig. 1 we show the phonon dispersion of MgB2, 
NbB2 and TaB2 for LA and optical E2g modes along 
T — A. For MgB2, a comparison of our calculations of 
E2g and LA phonon frequencies at T and A points with 
the previous calculations of Refs. |jl^, |l^, |l^ is given in 
Table I. As expected, our results are closer to the work 
of Ref. ||l5|. The difference in the i?2g phonon frequency 
at r, which has been found to be very sensitive to the 



Figure 2: The phonon linewidth 7qi, of MgB2 (upper panel), 
NbB2 (middle panel) and TaB2 (lower panel) for LA (solid 
bar) and optical i?2g (open bar) modes along F — A, calculated 
using the full-potential linear response method as described 
in the text. The phonon linewidth for the LA mode in MgB2 
has been multiplied by a factor of 10 for clarity. 



structural input and Brillouin zone integration jl^, arises 
due to the experimental lattice constants used in Refs. 
||l5| , |l6| , [l9| . The value calculated by Shukla et al. for 
E2g phonon frequency is somewhat lower. However, our 
calculated frequency for LA phonon mode at A is in good 
agreement with previous calculations. We expect ujqi, for 
NbB2 and TaB2 to have similar accuracy. 

In NbB2 and TaB2, the E2g phonon mode along T — 
A has considerably stiffened in comparison with MgB2. 
The E2g frequency in M gB2 changes from 78 to 72 meV, 
while for NbB2 and TaB2 it changes from 110 to 107 meV 
and 108 to 105 meV from T to A, respectively. The E2g 
frequency at F in TaB2, as calculated by Rosner et al. 
I^oj , is 98meV. The LA phonon mode at A for TaB2 
(22 meV) is almost half of that of MgB2 (40 meV) and 
NbB2 (37mey). Here, we Hke to point out that in the 
present work the number of q points chosen along T — A 
is not sufficient to say anything about the anomaly in the 
acoustical mode. 

The differences in the nature of electron-phonon inter- 
action between MgB2 and the transition-metal diborides 
NbB2 and TaB2 become quite apparent if one considers 
the electron-phonon contribution to the phonon lifetimes. 
In the case of MgB2, as shown by Shukla et al. 
the anharmonic effects Jit] , Q make negligible contribu- 
tion to the the phonon Hnewidth . Thus, the anomalous 
broadening of the E2g phonon linewidth along T — A un- 
derscores the strength of the electron-phonon coupling 
for this particular mode in MgB2 . To see what hap- 
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Figure 3: The partial electron-phonon coupling constant Aqi, 
of MgB2 (upper panel), NbB2 (middle panel) and TaB2 
(lower panel) for LA (solid bar) and optical i?2g (open bar) 
modes along F — ^4, calculated using the full-potential linear 
response method as described in the text. The total electron- 
phonon coupling constant Aq (lined bar) is also shown. 



Figure 4: The convergence of Uqv, ^y^v, and Aqi^ at A in 
the Brillouin zone for NbB^ as a function of k-points us- 
ing the double-grid technique as described in the text. In 
the present case the four sets correspond to 50 , 133, 133, and 
270 irreducible k-points for the electronic self-consistency and 
793, 793, 2413, and 5425 irreducible k-points for the Fermi 
surface sampling, respectively. 



pens in the transition-metal diborides, we show in Fig. 

2 the phonon linewidths of MgB2, NbB2 and TaB2 for 
LA and i?2g modes along T — A. For MgB2, our calcu- 
lated 7qi/'s are in reasonable agreement with the results 
of Ref. Note that the values shown in Fig. 2 cor- 
respond to twice the linewidth. From Fig. 2, it is clear 
that in MgB2 the electron-phonon coupling along T — A 
is dominated by the optical i?2g phonon mode with a 
maximum jE2g of 44 meV, and that the LA mode plays 
essentially no role. In contrast, in NbB2 and TaB2 (i) the 
linewidths are more than an order of magnitude smaller 
than in MgB2, for example the maximum ^E2g is only 
about 4 meV in TaB2 and (ii) the contribution from the 
i?2g mode decreases from 4 to 2.8 meV , while that due to 
LA mode increases from 0.02 to 1.5 meV, as one moves 
from r to A. The phonon linewidths of MgB2 and the 
transition-metal diborides NbB2 and TaB2, as described 
above, clearly demonstrate the differences in the strength 
and the nature of electron-phonon interaction in these 
systems. 

To see the strengths with which the LA and the i?2g 
phonon modes couple to the electrons, we show in Fig. 

3 the partial as well as the total electron-phonon cou- 
pling constant (Aq) along T ~ A for MgB2, NbB2 and 
TaB2. Certainly, the most striking feature of these sys- 
tems, as evidenced in Fig. 3, is the overall strength of 
the electron-phonon coupling in MgB2 (Aq ^ 7.9) as 
compared to NbB2 (Aq 1.5) and TaB2 (Aq ~ 2.5); 
nevertheless the additional feature of i?2g— dominated Aq 



in MgB2 giving way to L^-dominated Aq in NbB2 and 
TaB2 is just as striking. Thus, the electron-phonon cou- 
pling in the transition-metal diborides NbB2 and TaB2 is 
essentially due to the LA mode (the transverse acoustical 
mode makes some contribution), the contribution from 
the E2n mode being insignificant. For TaB2, Rosner et 
al. 12^ also found \E2g = 0-05 at F, in agreement with 
the present work. However, their ||2^ conclusion about 
the strength of the electron-phonon coupling in TaB2 is 
erroneous because it doesn't take into account the con- 
tributions from the acoustical modes properly. We also 
note that, in our opinion |pO|, [30|, the experimentally ||3^ 
deduced electron-phonon coupling in NbB2 and TaB2 
are underestimated. These differences in the electron- 
phonon coupling between MgB2 and the transition-metal 
diborides NbB2 and TaB2 may help explain why MgB2 
superconducts at 39 while NbB2 and TaB2 do not 
show any superconductivity down to 2i^. 

Finally, in Fig. 4 we show the convergence of Wq^, 
7qi/, and Xqi, at A in the Brillouin zone for NbB2 as a 
function of k-points using the double-grid technique as 
outlined in Ref. |Q. The use of double-grid technique 
allows one to construct two separate but commensurate 
k-grids, one for the electronic charge self-consistency and 
the other for Fermi-surface sampling. We employed 4 sets 
of double grids (i) (8, 8, 8, 24), (ii) (12, 12, 12 ,24), (iii) 
(12, 12, 12, 36), and (iv) (16, 16, 16, 48), where the first 
three numbers define the electronic self-consistency grid 
and the last number sets up the Fermi-surface sampling 
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rid which is commensurate with the first grid. We find [8] V. A. Gaspar ov et al, JETP Lett. 73, 532 (2001); cond- 



that i the results are converged for the (12.12.12.36) grid 



used in the present work. 

In conclusion, we have studied from first principles (i) 
the phonon dispersion, (ii) the phonon linewidth, and 
(iii) the partial electron-phonon coupling constant along 
T — A in MgB2, NbB2 and TaB2 in PG/mrnm crystal 
structure. We find that (i) in contrast to a strong and 
i?2£/-mode dominated electron-phonon coupHng in MgB2, 
the transition-metal diborides NhB2 and TaB2 have a 
relatively weak electron-phonon coupling which is dom- 
inated by the LA mode, (ii) the i?2g phonon linewidth 
is an order of magnitude larger in MgB2 than in NbB2 
or TaB2 , and (iii) the i?2g phonon frequency in NhB2 as 
well as TaB2 is considerably higher than in MgB2 while 
the LA phonon frequency at A for TaB2 is almost half 
of that of MgB2 or NbB2. 
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